Chaos and nights in the atom-photon interaction in cavity QED 
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We study dynamics of the atom-photon interaction in cavity quantum electrodynamics (QED), 
considering a cold two-level atom in a single-mode high-finesse standing-wave cavity as a nonlinear 
Hamiltonian system with three coupled degrees of freedom: translational, internal atomic, and the 
field. The system proves to have different types of motion including Levy flights and chaotic walkings 
of an atom in a cavity, ft is shown that the translational motion, related to the atom recoils, is 
governed by an equation of a parametric nonlinear pendulum with a frequency modulated by the 
Rabi oscillations. This type of dynamics is chaotic with some width of the stochastic layer that 
is estimated analytically. The width is fairly small for realistic values of the control parameters, 
the normalized detuning 8 and atomic recoil frequency a. ft is demonstrated how the atom-photon 
dynamics with a given value of a depends on the values of 8 and initial conditions. Two types of 
Levy flights, one corresponding to the ballistic motion of the atom and another one corresponding 
to small oscillations in a potential well, are found. These flights influence statistical properties of 
the atom-photon interaction such as distribution of Poincare recurrences and moments of the atom 
position x. The simulation shows different regimes of motion, from slightly abnormal diffusion with 
< x 2 >~ r 113 at 8 — 1.2 to a superdiffusion with < x 2 >~ r 2 " 2 at 6 = 1.92 that corresponds to 
a superballistic motion of the atom with an acceleration. The obtained results can be used to find 
new ways to manipulate atoms, to cool and trap them by adjusting the detuning 8. 

PACS numbers: 05.45.Mt, 42.50.Vk 



I. INTRODUCTION 

Cavity quantum electrodynamics is a rapidly develop- 
ing field of physics studying the interaction of atoms with 
photons in high-finesse cavities in a wide range of elec- 
tromagnetic spectrum, from microwaves to visible light, 
in such conditions under which both atoms and fields 
may manifest their quantum nature (for reviews on quan- 
tum and atom optics, see |IL 0). Modern experiments in 



cavity QED have achieved the exceptional circumstance 
of strong atom-field coupling for which the strength of 
the coupling exceeds both the atomic dipole decay and 
the cavity field decay providing manipulations with sin- 
gle atoms and photons (|, [| [|, 0, Trapped atoms 
and ions, interacting with laser fields in the regime of the 
strong coupling, have been used not only to study funda- 
mentals of quantum mechanics ||, ^| but also for appli- 
cations in the rapidly growing fields of quantum comput- 



2 



ing, quantum communication (see, for example, |Tc| , 
quantum chaos [jl2[ |l3| , and decoherence Jl4|, [ll| . 

A special comment should be related to the notion 
of chaos that we use in the paper. The real system is 
quantum and the quantum chaos per se does not exist. 
The notion "quantum chaos" just refers to the behavior 
of quantum systems whose classical counterparts behave 
chaotically. Through the paper the notion "chaos" will 
be used only in that meaning since all our consideration 
will be semiclassical. The key issue of quantum chaos is 
the correspondence between the classical and quantum 
pictures of chaotic dynamics. Chaotic dynamics in the 
atom-photon physics appeared in the paper p6| , where 
the semiclassical model of a single- mode, resonant, and 
homogeneously broadened laser, considered as an open 
dissipative system, has been shown to be equivalent to a 
Lorentz-type strange attractor, and the paper |l7|| , where 
Hamiltonian semiclassical chaos has been shown to arise 
in the Dicke model with non-resonant terms describing 
the interaction of an ensemble of identical two-level atoms 
with their own radiation field in an ideal resonant cavity. 
The fully quantum version of the latter model has been 
considered in It was shown in [18] that in a parame- 
ter range, corresponding to classical chaos, the evolution 
of the system becomes essentially quantum after the so- 
called "breaking time" Tn JlSj which obeys in this regime 
the logarithmic law r% ~ A -1 hi (const /Ti.), where A is the 
maximal Lyapunov exponent. All results presented here 
as chaotic dynamics are valid up to t < t%. 

The physical mechanism of semiclassical chaos in the 
Dicke model is tied to virtual transitions which are de- 
scribed by non-resonant (or counter-rotating) terms in 
the Dicke Hamiltonian. They are small under usual con- 
ditions. Trying to find another mechanism of local insta- 
bility in the Hamiltonian atom-photon dynamics with- 
out pump and losses, the authors of |2(], proposed 
the semiclassical model with atoms moving through a 
standing-wave cavity in a direction along which the cav- 
ity sustains a space-periodic field. The standing wave 



modulates the atom-field coupling providing in a certain 
range of the system's parameters intermittent Rabi oscil- 
lations. Another way to change the atom-field coupling 
is modulating the cavity length, studied in [^]. New ef- 
fects in the model with moving atoms may arise beyond 
the simple semiclassical approximation. Chaotic vacuum 
Rabi oscillations, a new kind of reversible spontaneous 
emission, have been shown p3| , pi| to occur in the model 
with interatomic quantum correlations. 

Experiments to study the quantum dynamics of clas- 
sically chaotic systems for the atom-photon interaction 
in cavities and traps have been intensively studied with 
cold atoms in a phase-modulated standing wave (2j| and 
in an amplitude-modulated standing wave |26|] (follow- 
ing to the proposition of (27)), and in a pulsed stand- 
ing wave [p8[ p9| . Cold sodium or cesium atoms, which 
are kicked by a periodically pulsed standing wave of far- 
detuned light, is an excellent experimental realization of 
a paradigm model of quantum chaos, a ^-kicked quan- 
tum rotor. Dynamical localization, that was observed in 
the atomic momentum distribution, is the quantum sup- 
pression of the classical momentum diffusion. A typical 
underlying phase-space structure of classical chaotic sys- 
tems consists of stochastic webs, islands, and chains of 
islands embedded in the stochastic sea. The chaotic mo- 
tion occupies a certain area in phase space. Because of 
islands and their boundaries, typical behavior of chaotic 
systems can be intermittent with long (quasi)regular os- 
cillations (the so-called Levy flights) interrupted by the 
chaotic pieces of trajectories. This intermittency leads 
to the anomalous diffusion with Levy distribution func- 
tions or a similar one, which have power-wise tails (for 
reviews on Levy processes in physics see |30|| , where the 
term "strange kinetics" was coined, and |3l|]). It was 
found with kicked cold cesium atoms j|2| that for certain 
pulses amplitudes, where the respective classical analog 
may exhibit anomalous diffusion, the momentum distri- 
butions were not exponentially localized for the time of 
observation (see also |}3L |34|, [3[| ) . It should be mentioned 
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that anomalous diffusion and Levy flights have been even 
earlier found with cold atoms and employed in a subrecoil 
laser cooling scheme |37], |3£| . 

In experiments with (5-kickcd atoms, the detuning be- 
tween the optical and atomic transition frequencies is 
large (relative to the natural line- width), so the prob- 
ability is small to find an atom, initially prepared in 
the ground state, in the excited state, and the excited 
state amplitude can be adiabatically eliminated |27|]. In 
this approximation, an effective Hamiltonian is that of 
a driven nonlinear oscillator with 3/2 degrees of free- 
dom. Generally speaking, the atom-photon interaction in 
a high-finesse cavity is, mainly, the interaction between 
internal (electronic) and external (motional) atomic de- 
grees of freedom and the cavity field. A corresponding 
one-dimensional model, including the interaction of all 
those degrees of freedom, has been introduced in papers 
p9| , piof in the context of Hamiltonian chaos. In this pa- 
per, we use a slightly generalized version of that model 
with three degrees of freedom to study the effects of chaos 
and Levy flights in the strongly-coupled atom-field sys- 
tem. 

In Sec. II the semiclassical basic equations in the 
Heisenberg representation are derived in the form of six 
coupled nonlinear equations with two control parameters, 
the normalized detuning 5 and atomic recoil frequency 
a. Taking into account that the frequency of the atomic 
(and field) Rabi oscillations is much more larger than the 
frequency of translational oscillations, reduced Bloch-like 
equations of the atom-field internal motion are derived 
and solved in Sec. III. The motion of the center of the 
atom mass is governed by a single equation of a paramet- 
rically perturbed pendulum. This dynamics generates a 
stochastic layer with an exponentially small width. In 
Sec. IV we simulate the basic set of equations at the 
fixed value of the normalized recoil frequency a = 10 
corresponding to a light atom in a microcavity with re- 
alistic parameters j§ H H [| • The maximal Lyapunov 
exponents and Poincare sections are calculated and it is 



found that the detuning S is a crucial parameter in tran- 
sition to chaos. Statistical properties of the atom-photon 
interaction are considered in Sec. V. We give an evidence 
of two types of Levy flights of an atom, one correspond- 
ing to an almost linear dependence of the atomic position 
x on time (superdiffusive and superballistic regime), and 
another one corresponding to small regular oscillations of 
the atom in the potential well. The Levy flights influence 
strongly such statistical properties of atoms and photons 
as distribution of Poincare recurrences and moments of 
x. The distribution of recurrences and time-evolution 
of the moments depend on the value of the detuning <5 
demonstrating different regimes from an almost normal 
diffusion to a superdiffusion. In Sec. VI we discuss briefly 
ways to manipulate the atomic motion by varying control 
parameters and initial internal atomic states. It is pos- 
sible, in particular, to cool and trap atoms by adjusting 
the detuning. 

II. BASIC EQUATIONS 

The basic model of interaction of radiation with mat- 
ter describes the energy exchange between a two-level 
atom and a single mode of the quantized radiation field 
in an ideal lossless cavity jy]. In general, this interaction 
should involve not only the internal atomic transitions 
and field states but also the center-of-mass motion of the 
atom. With the recoil effect to be included into consider- 
ation, the standard Jaynes-Cummings Hamiltonian can 
be extended as follows: 
P 2 

H = h fooJ a Sz + hu>fa + a — MIq{cl + +as + ) ■ cosktx , 

2m 

(1) 

where x and p are the atomic position and momen- 
tum operators, respectively. Transitions between two 
electronic states, separated by the energy huj a , are de- 
scribed by the spin operators with the commutation re- 
lations, [s + ,s_] = 2s z and [sz,s±] = ±s±. The pho- 
ton annihilation and creation operators with the com- 
mutation rule [a, d + ] = 1 characterize a selected mode 
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of the radiation field of the frequency Wf and the wave 
number kf in a lossless cavity of the Fabry-Perot type. 
The parameter fio is the amplitude value of the atom- 
field dipole coupling which depends on the position of 
an atom inside a cavity. As it is usually adopted in 
cavity QED, we write down the Heisenberg equations 
for the external atomic operators, x and p, and for 
slowly-varying amplitudes of the field and spin operators 
a(t) = aexp(— iui/t), a + (t) — a + exp(iu>ft), s±(i) = 
s±exp(±iujft) and s z (t) = s z 



d A - 
dt x ~ 


- £- 

TCI 1 


iP = 


= —hkftt (a + s_ + as + ) sin kfX , 




= —i(u>f — uj a )s + + 2iD,Qa + s z coskfX 


dt°- 


— i(u>f — u>a)s- — 2iQoa s z coskfX , 


dt ! ' z 


= — ii7 (a + s_ — as + ) coskfX , 




= — OIqs + cos k/x , 


Tt" = 


= iQ,Q§- coskfi . 



To avoid cumbersome notations in Eqs. (|2|) we use for 
the amplitudes the same notations as for the respective 
whole operators. 

In order to derive a tractable closed set of equations for 
expectation values from the Heisenberg operator equa- 
tions (Q), we use the semiclassical approximation. It 
means that all the operators and their products in Eqs. 
(^) are averaged over an initial quantum state, which is 
supposed to be a product state of the translational, elec- 
tronic, and field states. The expectation values of all the 
operator products are factorized to the products of the re- 
spective expectation values, e.g., ((a + s_ + as+) • sin ktx) 
= ((a + )(s_) + (a)(s+))sm(kf(x)). By choosing the fol- 
lowing dimcnsionless expectation values: x — kf(x), 
p = (p)/hk f , a x = i(a + a+), a y = ±(a - a + ), 
s x = + s+), s y = 2i(§- — s+), as dynamical vari- 

ables, we finally get from Eqs. (0) a nonlinear dynamical 



system 

x = ap , 

p = -2(a x s x + aySy) sin a; , 
s x = —6sy + 2a v s z cosx , 

(3) 

— & S ^iCI/rgSg COS iJj * 
Qi^g — Sy COS X ^ 
fty — COS X ■ 

that describes the interaction between three degrees of 
freedom in the strongly coupled atom-field system, trans- 
lational (p, x), electronic (s Xl s y ), and the field (a x ,a y ) 
ones. The dot in Eqs. (^) denotes the derivative with 
respect to the normalized time r = Clot- The con- 
trol parameters are the normalized recoil frequency a = 
hk 2 /mf2o and the normalized detuning between the fre- 
quencies of the field mode and the electronic transition, 
S = [ujj — uj a )/i} . The system (||) conserves the energy 

W = ^p 2 - 2(a x s x + aySy) cos x - Ss z , (4) 

and it possesses two additional first integrals 

4 + 4 + 4 = S 2 , 4 + a\ + s z = N . (5) 

The first one is simply the conservation of the atomic 
probability, and the second one is a conserved total num- 
ber of excitations, which is known to be a constant in 
the rotating-wave approximation. The equation of mo- 
tion for the atomic inversion s z is easily derived with the 
help of the integral N 

s z — 2(a x s y - a y s x ) cos x . (6) 

The semiclassical approximation, we used for (Q), means 
that the atom as a classical particle with external and in- 
ternal states moves in a self-consistent classical radiation 
field. 

From the dynamical systems point of view equations 
represent a system with 3 degrees of freedom (one 
degree of freedom per a canonically conjugate pair of the 
generalized momentum and coordinate) in 6-dimensional 
phase space. Indeed, we reinsert Eq.(^) for s z into (|^). 
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After that the system (||) describes 3 degrees of free- 
dom: the atom external coordinates (%,p), the atom 
internal coordinates (s x ,s y ), and the field coordinates 
(a x ,a y ). There are 2 constrains: energy integral E = 
H(p, X) s x , s y ; a x , a y ) and the spin integral S. The num- 
ber of excitations N should be used to determine s z as 
a function of other variables. That means that full dy- 
namics is defined in 4-dimcnsional hyperspace and should 
have domains of chaos due to its nonintegrability. One 
can say that the location of the domains of chaotic mo- 
tion, islands of regular dynamics, set of stationary points, 
and boundaries define a topology of the system's flow. 
The topology is two-parametric (E, S) and very compli- 
cated. Its description needs a separate investigation. 

III. REDUCED DYNAMICS AND THE 
ESTIMATION OF THE STOCHASTIC LAYER 
WIDTH 

We can simplify further the basic equations (|[) intro- 
ducing the combined atom-field variables 

u = 2(a x s x + ciySy), v = 2(a y s x - a x s y ) (7) 

and using the integrals (|J). As a result, one arrives at 
the closed five-dimensional dynamical system 

x = ap , 

p = —u sin a; , 

it = Sv , (8) 
v = -Su + 2{2Ns z - 3s 2 z + S 2 ) cos x , 
s z — —v cosx , 

which generalizes the corresponding equations of the pa- 
per jhj) (see Eqs. (3) therein which were derived in the 
limit of large N). It is obvious from Eqs. (||) that at 
exact resonance, 6 = 0, the slow translational variables 
x and p are separated from the fast atom-field u, v, and 
s z , and the system (||) becomes integrable. At S = 0, 
the atom moves in a spatially periodic optical potential 
U = — m(0)cosx with it(0) = u(t = 0) = const, and its 
center-of-mass motion satisfies the pendulum equation 



x + au(0) sinx = 0. It is easy to find that the dynamics 
of the internal atomic variable s z satisfies the following 
equation: 



s z = ±2 v / s z (s 2 z - Ns z - S 2 ) + C cosx 



(9) 



where C is an integration constant, and x(t) is a solu- 
tion of the pendulum equation mentioned above. The 
equation (^) can be integrated in terms of elliptic Ja- 
cobian functions with a solution which converges to the 
well-known Jaynes-Cummings semiclassical solution 
in the limit x = const. 

Out of resonance, at S ^ 0, the system (||) exhibits 
chaotic dynamics. In order to clarify the origin of chaos, 
consider Eqs. (^|) in the limit of large number of excita- 
tion TV and large detunings S comparing to S 2 . Taking 
into account that the normalized Rabi frequency is of the 
order of ~ \/N > 1 and is much more larger than the fre- 
quency of small translational oscillations, \fa <C 1, the 
equations for the fast atom-field oscillations are reduced 
to the Bloch-like form 



u = Sv , 

v = —Su + ANs z cos x 
s z = —v cosx , 



(10) 



where the function cos x may be considered approxi- 
mately as a constant c over a period of time of many 
Rabi oscillations. In this approximation, the quantity 



u A +v 2 + N(2s z y = R 



(11) 



plays a role of the length of Bloch vector, and the general 
solution of the Bloch-like equations (O) can be found 



u = u(0) 



N 



2c 



N 



s 



COS Q nT 



S 4NSc 
— v(0) sin Q n t + -tt~— s z (0)(1 - cos Q n t) , 



— — — u(0) sin SI a?t + v(0) cos JIatt- 
$2jv 



ANc 



N 



s z (0)smn N T , (12) 
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cS 

s z = u(0)-^2-(l — cos 51 att) — — — v(0) suiQnt- 

N 



C 

f2jv 



2c 



+ s z (0) 

where the quantity 

n N = ^5 2 + (2c) 2 N 



N I — — ] cos^att 

S2jV 



(13) 



is similar to the Rabi frequency. 

Since the function cos a; varies in time slowly compar- 
ing to the fast oscillating u, v, and z, the atom-field vari- 
able u may be considered approximately as a spatially 
independent frequency- and amplitude-modulated signal 
that parametrically excites the translational motion: 



x + au(r) sin x = , 



(14) 



that follows from the first two equations of the system 
(||). The modulation has especially simple form for ini- 
tial conditions u(Q) = v(0) = and s z (Q) =| S |, that 
corresponds to the atom prepared initially in the upper 
state while the field may be initially at any state, and 
c= 1: 



u(t) = (1 - cosier) 



(15) 



The Eq. (|14|) is derived from the following effective clas- 
sical Hamiltonian 
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H = —x — uj cosx + uo cosriArr-cosx = TLo + V , (16) 

where 7io is the unperturbed Hamiltonian of a free pen- 
dulum with the following frequency of small oscillations: 



u = —y/aN\6\ \ S\. 



(17) 



Rewriting the perturbation V = uj 2 cosUnt ■ cos a; in the 
form 

uj 2 

V = — [cos(x + 0.nt) + cos(x — Qnt)] , (18) 

one may consider ([l6]) as the Hamiltonian of a particle 
moving in the field of three plane waves in a frame moving 
with the phase velocity of the first wave, while the phase 



velocities of the second and third waves are Q,n and — Oat, 
respectively. 

As it follows from the general theory of perturbed mo- 
tion of Hamiltonian systems with 3/2 degrees of freedom 
|42| , the Hamiltonian ([l6|) induces chaotic dynamics in 
the so-called stochastic layer that appears due to the 
separatrix splitting. Let us consider the motion in the 
neighborhood of unperturbed separatrix of the pendu- 
lum (|l6|) . Consider the Poincare-Melnikov integral 



/oo 
{H ,V}dT 
-oo 



— uj 2 x s'm(x — Qnt)cIt , (19) 

J — oo 

where {H.q,V} is the Poisson bracket. This integral de- 
scribes changes of the atomic translational energy at the 
separatrix Ho = E s = uj 2 . To estimate ( |]~9|) for the dy- 
namics near the separatrix, one can use for x and x their 
known unperturbed separatrix solutions 



x s = 4 arctanexp[±w(r — t„ 
2oj 

x s = ±- 



(20) 



cosh[cj(r — t„)] ' 
where r„ is introduced as an initial condition. Using the 
solutions ([20]), we get 

1,00 dt ( n 



AE„ = ±2u z 



, Mil I r —t - 4>n) , (21) 

coshi \ uj ' 



where the new time t — lu(t — r„) and phase cj> n — £lNi~ n 
were introduced. The integral ( |2l"| ) has been calculated 
to give 



AE n = AE S sin 4> n , 



AE. = 2, n 2 e M^ N /2u) 



(22) 
(23) 



smh(iril?f /u>) 

The oscillating function AE n has simple zeroes that im- 
plies transversal intersections of stable and unstable hct- 
eroclinic manifolds of saddle points known as a com- 
plicated heteroclinic structure. On the basis of general 
properties of motion near the separatrix, the separatrix 
map can be introduced [E2| . 



E n+1 = E n + AE s sm<j> n , (24) 
t>n+l = 4>n H hi -r= . (25) 
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The condition 



K = max 



exp(7rilAr/2o;)| sin <p n \ 



> 1 (26) 



uj \E n+ i — E s \ smh(irQ,N /u) 
defines the stochastic layer width which can be estimated 
in the case of the large parameter Un/uj 1 (see the 
respective estimations with real atoms in the concluding 
section) as follows: 



SE S 



\E, 



n+l 



exp 



(27) 



uj \ 2uj 

The dimensionless width of the stochastic layer is finally 
given by 

3 

exp 



5E S 

>2 



8tt 



n 



N 



u>* \ UJ 

where the large parameter 

n N S 2 + AN 



N 

2uj 



(28) 



w 2y/aN | S 

under the conditions N,S S \ 
estimated as 



S 



(29) 



UJ 



N 



1 and S 2 < AN is 



(30) 



For the considered case the width of the stochastic layer 
of the reduced atom- field dynamics is exponentially small 
in 



), multiplied by a large parameter. Due to (30) the 
final width depends on the control parameters N, a and 
S, and the formula ((2^) is useful in estimating the ranges 
of the control parameters where one may expect chaotic 
motion. 

The estimation ( p8| ) provides the lower bound for the 
the width of the stochastic layer that appears due to the 
simplest harmonic modulation ( |l5| ) . Small changes in en- 
ergy produce comparatively small changes in frequency 
of oscillations. Nearby the bottom of potential wells and 
high over potential hills (where the energy is much less 
and much greater than E s ), small changes in frequency 
give rise to respectively small changes in phase during the 
period of oscillations. Nearby the unperturbed separa- 
trix, where the period goes to infinity, even small changes 



in frequency lead to dramatic changes in phase. This is 
the reason of exponential instability of the parametric 
oscillator ( |l4| ) and (|l5|) which models chaotic motion of 
the atom moving through a periodic standing wave. 



IV. LYAPUNOV EXPONENTS AND POINCARE 
SECTIONS 

In this section, we present numerical simulations with 
the basic set of equations (Q) and the integrals of motion 
(|) and (|) with S 2 = 3/4 (the actual value of S 2 has no 
importance since it always can be renormalized to 1) and 
N = 10. The system @ has two control parameters, the 
normalized detuning S between the atomic transition and 
cavity frequencies and the normalized recoil frequency a. 
As it will be estimated in the concluding section, a is in 
the range from 10~ 5 to 10 -2 for real atoms. We choose 
a = 10~ 3 in simulations throughout the paper. 

The detuning 5, as it was shown in m 0, is the 
crucial parameter in transition to chaos in the atom-field 
system with the center-of-mass motion. It is obvious from 
the set (g), which is equivalent to the basic one (§), that 
at exact resonance, S = 0, the motion is regular. At 
large detunings, 6 3> N, the motion is expected to be 
quasiregular since the nonlinear term in the fourth equa- 
tion of the set (||) is small, compared to the linear term 
of the same equation. With far-detuned light, one does 
not expect pronounced atomic Rabi oscillations. In order 
to find the range of the detunings, where the motion is 
expected to be chaotic, we compute the dependence of 
the maximal Lyapunov exponent A on the detuning S. 

Lyapunov exponents characterize the behavior of close 
trajectories in phase space. Consider a trajectory, some 
sequence of time instants tq, t\, t%, ... with equal in- 
tervals At, and a ball around the initial point of the 
trajectory. Lyapunov numbers Aj (J = l,...,M; M is a 
number of variables) show per-interval At changes of the 
axes of the "ellipsoid" of the deformed ball [[l6) . In our 
case M = 6 after exclusion of s z . The A:-th Lyapunov 
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FIG. 1: The maximal Lyapunov exponent A in units of the 
maximal atom-field coupling rate Qo versus the atom-field 
detuning S in units fl . (a) s z (0) = -0.863. (b) s z (0) = 0. 
(c) s z (0) = 0.863. 

exponent is defined as Afc = In Afc. Typically Afc depends 
on time and Afc should be replaced by their mean values 
p6[ . In Hamiltonian systems, due to the phase volume 
conservation, Ai • A2 • ... ■ Am — 1 and Ai + ... + Am = 0. 
For integrable system all Afc are pure imaginary and they 
make pairs: A2 = — Ai, A4 = —A3, Ag = — A5, since the 
number M of equations is even. This result follows from 
the so-called Liouville- Arnold theorem ||47fl. In our case, 



for the reduced system of 6 variables (p, x; a x , a y \ s x , s y ) 
and two constrains (integrals of motion), we have two 
imaginary pairs, say A1.2 = ±«cri, A3,4 = ±z<72 (ci,2 real), 
and A5, A6 that satisfy the condition A5 + Ag = 0, and 
therefore A5 • Ag = 1, i.e. Ag = I/A5. The chaos means 

If, say A 5 < 1 (A 5 < 0), then 



that A5 g are real 
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Ag > 1 (Ag > 0) and Ag is called maximal Lyapunov 
exponent. It has a nice physical meaning: the maximal 
Lyapunov exponent measures a rate of the separation of 
initially close trajectories, and typically for the practi- 
cal goal, the mean value A = Ag over time is used. To 
compute A, we use the standard algorithm M] 



A 



lim — In -— - k \ , 

„_>no n ^ — ' A(Tfc_l) 



(31) 



fc=l 



where A(rfc_i) is a distance between two close trajecto- 
ries at time Tfc_i, and the value A(rfc)/A(rfc_i) shows the 
level of separations of the trajectories during the interval 

(Tfc_l,Tfc). 

In the case that separation doesn't go exponentially, 
A = 0. This happens when 5 = since the system be- 
comes integrable, and the exponential separation disap- 
pears. 

The corresponding results of computing maximal Lya- 
punov exponent A are presented in FIG. |l| with three 
different initial values of the atomic population inver- 
sion s z (0) = —0.863; 0; 0.863, respectively. The other 
initial conditions are the following: a x (0) = s x (0) = 0, 
x(0) = 3.14, p(0) = 2, and a y (0) and s y (0) are found 
from the Eqs.(|4|) and (||) with given S = and 
N = 10. The value s z (Q) = —0.863 corresponds to the 
atom initially prepared closely to its ground state for 
which s z = — y/3/2. Note that the unusual amplitude 
values of the atomic population inversion we have are 
the result of the chosen normalization S 2 = 3/4. The 
atom with s z (Q) = 0.863 is prepared closely to its ex- 
cited state. In both the cases, the initial components of 
the transition electric dipole moment, s x (0) and s y (0), 
are almost zero with |s z (0)| = 0.863. The atom with 
s z (0) = has a maximal electric dipole moment. 
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FIG. 2: Projections of the Poincare sections on the plane of the atomic momentum p in units hkf and the position in units of 
kj 1 . (a) s x (0) = -0.8660254 and 5 = 1.2. (b) s z (0) = -0.863 and 8 = 1.2. (c) s z = -0.8660254, 8 = 1.92, and p(0) = 2.1 in 
all the above fragments, (d) s z = —0.863, 5 = 1.92, and p(0) = 20. x,p are dimensionless. 



As it was expected, at exact resonance (8 — 0), the 
maximal Lyapunov exponent is exactly equal zero in all 
the cases assuming a regular motion. As it follows from 
the results of previous section, a stochastic layer ap- 
pears with infinitesimally small values of detuning, but 
its width decreases fast with increasing 5 (see (|28|) and 
(p9|)). We find A ~ with |<5| > 3. In physical terms, 
it means that at exact resonance an atom will periodi- 
cally exchange excitation with the field, whereas far off 
resonance its internal states will not (almost) be affected 
by the field. This interplay results in a maximum of the 
A (8) dependence with almost the same maximal values 
for s z (0) — ±0.863 and 0. The results, however, are 
different in the range —3 < 8 < 3 with different initial 



values of the population inversion s z . It is easy to under- 
stand why it is. As it follows from the Bloch-like solution 
( |l2| ) for s Z) the atom starting, say, in its ground state, 
s 2 (0) = — v3/2 and u(0) = v(0) = 0, could reach the up- 
per state s z = y/3/2 only with 6 = 0. The same is valid 
with the other initial values of s z (0): the atom starting 
with s z (0) = \/3/2 (or with s 2 (0) = 0) will not reach 
s z = — V3/2 (or s z = ±V3/2) except for the case of ex- 
act resonance, 8 — 0. Thus, the dependencies X(S) are 
different in the range — 3 < 8 < 3 with different values 
of s 2 (0) in spite of the fact that the maximal Lyapunov 
exponent is computed over the rather long trajectory. 

The model Hamiltonian (0) can be easily generalized to 
an ensemble of indistinguishable two-level atoms. In the 
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semiclassical approximation, we have not observed any 
pronounced differences in the strength of chaos (that is 
characterized by the values of A) with different initial in- 
ternal atomic states. Interatomic quantum correlations, 
which occur through the mediation of the field gener- 
ated by the atomic ensemble, have been shown in H] 
(where the model with hot moving atoms but without 
recoil has been considered) to play a significant role in 
the atom-field dynamics. Much more strong chaos has 
been numerically found [p4j in the vacuum Rabi oscilla- 
tions with atoms initially prepared in the so-called super- 
fluorescent state (with all the atoms to be uncorrelated 
initially and occupying their excited states) than with 
atoms initially prepared in the superradiant state (with 
initially strongly correlated atoms having a macroscopic 
electric dipole moment). 



We numerically construct single-trajectory Poincare 
sections of motion in the system (|^) with three degrees 
of freedom and project them on the plane of the atomic 
external variables (x,p). FIG. || presents these sections 
with the atom initially prepared close to the ground state 
and with two different values of the atom-field detuning 
S = 1.2 (see (a) and (b) fragments) and 8 — 1.92 (see 
(c) and (d) fragments). In the latter case chaos is not as 
strong as in the first case (see FIG. |l|a). A fairly regular 
web structure, that is seen in the fragment (b) computed 
over a comparatively short integration time, breaks down 
with increasing integration time (the fragment (a)). For 
comparison, we present in FIG. ^|c the Poincare section 
computed under the same conditions as in FIG. |^a but 
with S = 1.92 (an additional trajectory with p(0) = 0.2 
is plotted in the fragment (c)). FIG. ||d demonstrates 
the Poincare section at 6 = 1.92 with an increased initial 
momentum p(0) = 20. 



a 




X 



FIG. 3: The same as in FIG. | with s z (Q) = -0.863 and 
p(0) = 1. Change of density appears after (a) At = 75, 000 
and S = 1.2 and (b) At = 49,000 and S = 1.92. x,p are 
dimensionless. 

V. STATISTICAL PROPERTIES OF THE 
ATOM-PHOTON INTERACTION 

In this section a detailed analysis of the atom and pho- 
ton chaotic dynamics will be considered. A sensitive con- 
trol parameter is S = (ujf — w Q )/f2o, detuning of the field 
and atom frequencies. For the sake of convenience we 
specify two values of 5: 1.2 and 1.92. It follows from cal- 
culating maximal Lyapunov exponents in the previous 
section that the smaller is S (in the range 0.5 <| 5 |< 3), 
the stronger is mixing and chaos, and one can expect 
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FIG. 4: Levy flights of an atom in a cavity, (a) A long ballistic 
flight at S = 1.2 is evident, (b) Ballistic flights intermittent 
with stagnation phases of motion at S = 1.92 are shown, (c) 
Two different types of the Levy flights in the plane x — p. 
Time is in units of fig -1 . x,p, r are dimensionless. 



that the case with 5 — 1.92 is, being chaotic (but not 
with the atom prepared in the upper state or close to 
it), more intermittent than the case 6 = 1.2. This prop- 
erty of the atom-photon dynamics will be quantitatively 
characterized below. 

The difference of a trajectory projection on the x — p 
plane [a x = s x = 0) is evident from FIG. |^ where the den- 
sity modulation has been used: a change of each density 
appears after At points of the mapping the trajectories 
(At = 75,000 for (a) and At = 49,000 for (b)). The 
narrow strips of the same density in FIG. [|b indicate a 
long stay of atom in the corresponding part of the x — p 
plane with oscillations in the potential well and a small 
change of the amplitude of the oscillations. In contrast 
to this pattern, the distribution of densities in FIG. ||a is 
more uniform manifesting much better mixing, although 
some traces of the intermittency persist. 

The difference between 5 = 1.2 and 5 — 1.92 is also 
evident from FIG. ^ where a dependence x — x(t) is 
shown. The intermittent case (S = 1.92) has very long 
"flights" known also as Levy flights |30|, [n]. There are 
two types of flights in FIG. ||b. One category of flights 
corresponds to the almost linear dependence of x = x(r), 
while the other one corresponds to the stagnation of the 
trajectory near some value of x. FIG. |^c shows the flight 
in the x—p plane where the ballistic dynamics coexists or 
alternate stagnations. Both categories of flights are well 
understood from FIG. §c,d: ballistic dynamics along x 
in FIG. ||c is responsible for the linear dependence of 
x = x(t), while the trajectory can stay very long near 
the saddle points as in FIG. §1 (the dark area near a 
saddle point). The case of 8 = 1.2 in FIG. ||a,b is very 
different and flights of both categories are rare, if ever. 

Just the presence of flights and intermittent behavior 
of the physical variables strong fluences the statistical 
properties of atoms and photons. We will use two impor- 
tant characteristics of the atom-photon variables: distri- 
bution of Poincare recurrences and moments of the atom 
coordinate x. Consider a small phase volume AT and 
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P(AT; t) as a probability density of a trajectory to re- 
turn first time back to Ar at time instant r £ (r + dr) 
if initially started at Ar at r = 0. Then the density 
probability to return first time to Ar is 



P(t) = lim — P(AT; r) 
v ' Ar^o Ar v ; 



(32) 



with a normalization condition 

/ P(r)dT = lim — / P(AT; 
Jo V ' Ar^o Ar J Q V 



T)dT = 1 . (33) 



The probability P{t) does not depend on the choice of 
Ar and for "good" chaotic mixing decays exponentially 



P(t) = (l/h)e~ hr 
with the mean recurrence time 



POO 

l/h= / TP(r)dT 
Jo 



(34) 



(35) 



and h as Kolmogorov-Sinai entropy. 

The general situation is more complicated since an al- 
gebraic behavior 



P(t) ~ 1/r 7 



(36) 



is possible for large r due to intermittent chaos, ror 
bounded Hamiltonian dynamics, r rcc < oo (Kac lemma), 
and the condition 7 > 2 should exist. Nevertheless, 
strongly intermittent dynamics sometimes docs not per- 
mit us to achieve the limit at r — > 00 and many different 
intermediate asymptotics can appear. HG. || shows the 
distribution of recurrences that is close to the exponen- 
tial one as in (p4[) for 5 = 1.2 and to the algebraic one as 



in Q with 7 > 2 (t > 10 5 ) for S = 1.92. 

The difference due to intermittency also occurs for the 
moments 



„2m\ 



u(m) 



(37) 



where the so-called transport exponent fi(m) varies for 
different time windows. The behavior of (x 2m ) is shown 
in FIG. |. For 5 = 1.2, the value fi(l) is close to 2 
for r < 10 3 and corresponds to the ballistic dynamics. 





FIG. 5: Time evolution of the 2mth-order moments of the 
atom position a; on a logarithmic scale with the values of 
slopes indicated for each m. The insets show the respective 
distributions of the Poincare recurrences, (a) S — 1.2. (b) 
8 — 1.92. x,t axe dimensionless. 



For r > 10 3 and 6 = 1.2, we have ^(1) ps 1.13 that 
corresponds to a weak superdiffusion that is fairly close 
to the normal diffusion with jtx(l) = 1. A very different 
behavior for moments appears for 6 = 1.92 where there 
are many long-lasting flights. For t > 10 3 , fi(l) ps 2.2 
that corresponds to a superballistic transport with an 
acceleration. This behavior can be explained as a result 
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of long flights when atoms move in the photon's field 
acquiring acceleration. This type of transport is self- 
similar and fi(4) sa 8.8 = 4/i(l). 

VI. MANIPULATION OF ATOMS 

In this section we would like to make a few comments 
related to the manipulation of atoms by changing differ- 
ent control parameters. As it was shown in Sec. V, a 
change of 6 leads to a possibility of a sensitive control of 
the Levy flights and, as a result, to cool the atoms which 
have the lower chaotic dispersion the longer the flight is. 
At the same time, simulations show fast mixing on the 
&X and s x — s y planes. More precisely, spectral prop- 
erties of the atomic dynamics are sensitively controlled 
by the parameter 8. Let us demonstrate it using a sim- 
plified analysis. 

Consider x = x(t) as the only variable that describes 
the dynamics or the most essential part of the dynamics, 
and introduce a generation function 

G{x, t; v) = e *"[*M-»o] j XQ = 2,(0) _ ( 38 ) 

Then 



I(t;x ) 



dv G(x, t; v) = 2it5[x{t) ~ x ] . (39) 



The expression I(t; xq) can be "coarse-grained" over xq, 
i.e. 



(I(t;x )) Xo 



Ax 



dx J(r, x ) 



2tt 
Ax 



dx 8[x(t) - Xq] . (40) 



The presence of 5-function indicates recurrences to xq 
within an interval Axq at time instant r within an inter- 
val Atq. For r — > oo we can neglect n-triple recurrences 
from n > 2 and leave only the first recurrences. Then 



(I(t;x )) X0 =P(t) 



E 

XqGAxq 



const / 



dx(r) 



dxo 



(41) 



The expression ([IT]) shows that for "good" chaotic sys- 
tems 



\dx(r) / dxo \ ~ expQir) 



(42) 



and we arrive at (pj). For the intermittent dynamics, the 



sum in (41) consists of two types of terms, the same as 
and the algebraic growth 



\dx/dxo\ 



(43) 



with the value of 7 depending on the type of intermit- 
tency. For a fairly large r, the term (^) survives and we 
arrive at the case (p6[). 
From another side, 

Q2m 



■(G(x, r; v) v=a = (~l) m (\x(r) - x \ 2m ) , (44) 



and we obtain the moments of x(r). This shows that the 
moments and their spectral properties are coupled to the 
recurrences distribution through the generating function 
G(x, r; v) which one would expect to obtain from exper- 
iments. When the moments are infinite, the expression 
(H) can be replaced by the following one 



(G{af,T; u p ) v=0 = const (\x(t) - x f) , (45) 



x(t)=xo 



with an appropriate value of (3 (see more discussion in 
II)- 

The main way of controlling the properties of G(x, t; v) 
is to change the system's topology in phase space. Speak- 
ing about the topology, we have in mind the phase 
pattern (see also the end of the Sec. II) that includes 
the singular points, curves, and partitioning of the do- 
mains of chaos and islands. To illustrate how the sys- 
tem is sensitive to small variations of the initial condi- 
tions that change the full energy, we show in FIGs. ||a 
and b the Poincare sections with s z (0) = 0.863 and 
s 2 (0) = 0.8660254, respectively, at 6 = 0.4 and under the 
other equal conditions. Very small difference in the values 
of the initial tipping angle between the direction of the 
Bloch vector and the axis z gives rise to cardinally differ- 
ent motion with p(0) = 2, chaotic oscillations in the wide 
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FIG. 6: The same as in FIG. | but with p(0) = 2, 5 = 0.4 
and the atom initially prepared nearly in its excited state 
with two slightly different values of the population inversion, 
(a) s z (0) = 0.863. (b) s*(0) = 0.8660254 (for comparison, 
additional trajectories with p(0) > 2 are shown). x,p are 
dimensionless. 



range of the atomic momenta with s z (0) = 0.863 and 
small regular translational oscillations nearly the bottom 
of a potential well with s z (0) = 0.8660254. Transition 
from order to chaos takes place with the latter value of 
the initial atomic population inversion only with much 
more large values of the initial momentum, p(0) > 40. 

The original system (|^) has three degrees of freedom 




5 = 0.225 




1 57 



471 6.2 1.57 



FIG. 7: Bifurcation of the saddle point into the elliptic point 
on the plane x — p under changing the detuning S. x,p are 
dimensionless. 

and it is not studied yet. Nevertheless, we were able to 
demonstrate by simulation a bifurcation of a hyperbolic 
point into the elliptic one, although we are not able to 
provide an analytical description at the moment since 
the FIG. ^ is just a projection of a trajectory in the 
4-dimensional hyperspace onto the plane (p, x) . By a 
change of 8 near 8* ~ 0.222, the saddle on the x — p plane 
with a x = s x = transforms into the elliptic point. The 
trapping potential well of the finite size on x — p plane 
occurs as a result of the bifurcation: this bifurcation will 
be studied in detail in another paper. 

VII. CONCLUSION 

A system with one or more cold atoms strongly cou- 
pled to a single mode of the cavity field is ideal for testing 
fundamentals of quantum mechanics and its correspond- 
ing to classical mechanics. Based on our understanding 
of the nonlinear dynamics of the atom-photon interaction 
in a standing-wave high-finesse cavity, new ways to ma- 
nipulate and control atomic motion can be opened. We 
have shown that the motion is very sensitive to the atom- 
cavity detuning 8. Varying 8, one can design topology 
of the underlying phase space creating zones of trapping, 
quasi-trapping or acceleration, quasi-regular and stochas- 
tic webs, etc. It may provide new schemes for cooling, 
trapping, and accelerating atoms. 

To give an idea about the values of the magnitudes we 
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have used in numerical simulations, we need to estimate 
the range of values of the normalized recoil frequency 
a = hk^/mflo with real atoms and cavities. We will use 
the parameters of the real experiments with single atoms 
in the strong-coupling regime |3[ for which the max- 
imal atom-field coupling strength Hq exceeds the decay 
rates of the cavity field and of the atomic dipole. Atoms 
were collected in a magneto-optical trap and cooled down 
to microKelvin temperatures, before entering a micro- 
scopic high-finesse Fabry-Perot cavity with Q ~ 10 6 , 
fl ~ 2tt(10 7 - 10 8 ) Hz and k f ~ 2tt • lO 6 ™" 1 . With 
these values of the parameters, one can estimate a to be 
in the range 1CP 5 — 1CP 2 depending on the atomic mass 
and fig. 
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